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Abstract 

We study the relativistic electron-positron field at positive tem- 
perature in the Hartree-Fock-approximation. We consider both the 
case with and without exchange term, and investigate the existence 
and properties of minimizers. Our approach is non-perturbative in 
the sense that the relevant electron subspace is determined in a self- 
consistent way. The present work is an extension of previous work by 
Hainzl, Lewin, Sere, and Solovej where the case of zero temperature 
was considered. 
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Introduction 

In Coulomb gauge and when photons are neglected, the Hamiltonian of 
Quantum Electrodynamics (QED) reads formally [2J [T^l [T5] 

W = j tt*(z)Z>°tf (a;) dx- j <p(x)p(x)dx + ^JJ n~T^ dx d V • C 1 ) 

Here ^(x) is the second-quantized field operator satisfying the usual anti- 
commutation relations, and p(x) is the density operator 

P[ ' 2 2 

(2) 

where cr is the spin variable. In (pQ), D° = — ia • V + is the usual free 
Dirac operator, a is the bare Sommerfeld fine structure constant and <p> is 
the external potential. The matrices a = (01,02,03) and (3 are the usual 
4x4 anti-commuting Dirac matrices. We have chosen a system of units 
such that H = c = m = 1. In QED, one main issue is the minimization 
of the Hamiltonian (pQ). However, even if we implement an UV-cutoff, the 
Hamiltonian is unbounded from below, since the particle number can be 
arbitrary. 

In a formal sense this problem was first overcome by Dirac, who sug- 
gested that the vacuum is filled with infinitely many particles occupying 
the negative energy states of the free Dirac operator D . With this axiom 
Dirac was able to conjecture the existence of holes in the Dirac sea which 
he interpreted as anti- electrons or positrons. His prediction was verified by 
Anderson in 1932. Dirac also predicted [6, 7J the phenomenon of vacuum 
polarization: in the presence of an electric field, the virtual electrons are 
displaced and the vacuum acquires a non-uniform charge density. 

In Quantum Electrodynamics Dirac's assumption is sometimes imple- 
mented via normal ordering which essentially consists of subtracting the 
kinetic energy of the negative free Dirac sea, in such a way that the kinetic 
energy of electrons as well as positrons (holes) becomes positive. With this 
procedure the distinction between electrons and positrons is put in by hand. 
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It was pointed out in [H] (see also the review [13]), however, that normal 
ordering is probably not well suited to the case a 7^ of interacting particles 
(the interaction is the last term of ([T])). Instead a procedure was presented 
where the distinction between electrons and positrons is not an input but 
rather a consequence of the theory. The approach of [H] is rigorous and fully 
non-perturbative, but so far it was only applied to the mean-field (Hartree- 
Fock) approximation, with the photon field neglected. It allowed to justify 
the use of the Bogoliubov-Dirac-Fock model (BDF) [3], studied previously 
in [TO] QTJ [p2]. The purpose of the present paper is to extend these results 
to the nonzero temperature case. 

The methodology of [13] is a two steps procedure. First, the free vacuum 
is constructed by minimizing the Hamiltonian ([T]) over Hartree-Fock states in 
a box with an ultraviolet cut-off, and then taking the thermodynamic limit 
when the size of the box goes to infinity. The limit is a Hartree-Fock state 
describing the (Hartree-Fock) free vacuum |14(, [13] . It has an infinite energy, 
since it contains infinitely many virtual particles forming the (self-consistent) 
Dirac sea. We remark that this state is not the usual sea of negative electrons 
of the free Dirac operator because all interactions between particles are taken 
into account, but it corresponds to filling negative energies of an effective 
mean- field translation invariant operator. 

The second step of [H] consists of constructing an energy functional that 
is bounded from below in the presence of an external field, by subtracting the 
(infinite) energy of the free self-consistent Dirac sea. The key observation 
is that the difference of the energy of a general state P minus the (infinite) 
energy of the free vacuum V°_ can be represented by an effective functional 
(called Bogoliubov-Dirac-Fock (BDF) [3J) which only depends on Q = P — 
V°_, describing the variations with respect to the free Dirac sea. The BDF 
energy was studied in |101 \TT\ [12] . The existence of ground states was shown 
for the vacuum case in [1U|, [TT] and in charge sectors in [p2] . For a detailed 
review of all these results, we refer to [13] . An associated time-dependent 
evolution equation, which is in the spirit of Dirac's original paper [6], was 
studied in [15] . 

Let us now turn to the case of a non zero temperature T = 1//3 > 0. We 
consider a Hartree-Fock state with one-particle density matrix < P < 1. 
Because of the definition of the Hamiltonian ([1]) and the anticommutator in 
([2]), it is more convenient to consider as variable the renormalized density 
matrix 7 = P — 1/2. We remark that the anticommutator in ([2]) is a kind of 
renormalization which does not depend on any reference as normal ordering 
does (it just corresponds to subtracting the identity divided by 2). The 
anticommutator of ([2]) is due to Heisenberg [16] (see also [181 Eq. (96)]) and 
it is necessary for a covariant formulation of QED, see |22} Eq. (1.14)] and 
m Eq- (38)]. 

Computing the free energy of our Hartree-Fock state using ([!]) (and 
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ignoring infinite constant terms) one arrives at the following free energy 
functional O [3] 



^ EC ( 7 ) = trp» 7 )- a /^K(x ) + f // b^aM 

a f f tr C 4 \j(x,y)\ [ 
2 J J \x-y\ 

where the entropy is given by the formula 



rs( 7 ) (3) 



S( 7 ) = - tr (( J + 7 ) ln(| + 7 )) - tr ((§ - 7) M§ - 7)) • (4) 

The (matrix- valued) function j(x,y) is the formal integral kernel of the 
operator 7 and p-y(x) := ti£4 j(x, x) is the associated charge density. The 
above formulas are purely formal; they only make sense in a finite box with 
an ultraviolet cut-off, in general. 

As in [H] the first step is to define the free vacuum at temperature T, 
which is the formal minimizer of © when <p = 0. Following [14] , one can 
first confine the system to a box, then study the limit as the size of the box 
goes to infinity and identify the free vacuum as the limit of the sequence of 
ground states. Alternatively, it was proved in [14J that the free vacuum can 
also be obtained as the unique minimizer of the free energy per unit volume. 
In the nonzero temperature case, this energy reads 

T T (j) = 

1 1 iwPiFWj.ift-A// trc : l7W 7, (9)1 ^„ 



(27r) 3 Jb(0,A) (2T) 5 iii!(0,A)2 \P- 

+ 7^3 / M (I + 7 (p)) In (i + 7 (p)) + (I - 7 (p)) In (i - 7 (P))] dp 
K Z7r ) Jb(o,a) 

and it is defined for translation-invariant states 7 = j(p) only, under the 
constraint —1/2 < 7 < 1/2. Here, B(0, A) denotes the ball of radius A 
centered at the origin. The real number A > is the ultraviolet cut-off. We 
shall prove in Theorem H] that the above energy has a unique minimizer 7 , 
and prove several interesting properties of it. In particular, we shall see that 
it satisfies a nonlinear equation of the form 

^ = 2 (i + e ^o " i + e -/"v) (5) 

i.e. it is the Fermi-Dirac distribution of a (self-consistent) free Dirac opera- 
tor, defined as 

D, =D°-ap*4. 

\x — y\ 

(The last term stands for the operator having this integral kernel.) This 
extends results of HH to the T > case. 
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The next step is to formally subtract the (infinite) energy of 7 from the 
energy of any state 7. In this way one obtains a Bogoliubov-Dirac-Fock free 
energy at temperature T = 1/(3 which can be formally written as 

^• r ( 7 ) = "^Q ED ( 7 )-^Q ED (^)" 

= T% 7 u )-a J (p{x)p [7 _^ 0] {x) + - J J j— — 

g f f ^\(l-l°)(x,y)\ 2 (6) 
2 J J \x — y\ 

where H is the relative entropy formally defined as 

Ttf (7,7°) = "tr(^ ( 7 - f)) - TS( 7 ) + TStf°)\ 

We shall consider external field of the form ip = v * A, where v represents 
the density distribution of the external particles, like nuclei, or molecules. 

In Section [2.21 we show how to give a correct mathematical meaning to 
the previous formulas and we prove that the BDF free energy is bounded 
from below. An important tool is the following inequality 

TF(7,7°)>tr [1^,1(7 -7 ) 2 ] > tr [|Z?°| (7 - 7°) 2 ] • (7) 

This implies that the relative entropy can control the exchange term and 
enables us to show that Tt is bounded from below. 

Unfortunately, like for the T = case, the free BDF energy is not con- 
vex, which makes it a difficult task to prove the existence of a minimizer. 
Although we leave this question open, we derive some properties for a po- 
tential minimizer in Section [2.2i In particular we prove that any minimizer 
7 satisfies the following nonlinear equation 

7 = 2 (l + e^-r ~ 1 + e-^-r ) (8) 
where the (self-consistent) Dirac operator reads 



_D 7 = D° + a/) 7 * I • I 1 — a<f — a- 



x-y\ 



Compared with the zero temperature case, the main difficulty in proving 
the existence of a minimizer comes from localization issues of the relative 
entropy which are more involved than in the zero temperature case. 

As a slight simplification, we thoroughly study the reduced Hartree Fock 
case for T > 0, where the exchange term (the first term of the second line 
of ©) is neglected. In the zero-temperature case, this model was already 
studied in detail in [11] and [Sj. The corresponding free vacuum is now 
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simple: it is the Fermi-Dirac distribution corresponding to the usual free 
Dirac operator D°, 



7 



2ll + e^ D ° l + e-P D ° 



The reduced Bogoliubov-Dirac-Fock free energy is obtained in the same way 
as before by subtracting the infinite energy of the free Dirac see 7 to the 
(reduced) Hartree-Fock energy. It is given by 



^? d ( 7 ) = TF( 7 , 7 °) -a J W[7 _ 7 o] + f 



\x - y\ 



H (7, 7 ) being defined similarly as before. As this functional is now convex, 
we can prove in Theorem [2] that it has a unique minimizer 7, which satisfies 
the self-consistent equation 



1/1 1 

7 



2 V 1 + eP Ds > 1 + e-P D ~-i 
where 

:= D + apj_ y o * \ ■ I" 1 — atp 

in this case. 

Additionally we show in Theorem [3] that this minimizer has two interest- 
ing properties. First, 7— 7 is a trace-class operator. In the zero temperature 
case, on the other hand, it was proved in [11] that the minimizer is never 
trace-class for a > 0. This was indeed the source of complications concern- 
ing the definition of the trace (and hence of the charge) of Hartree-Fock 
states |10| when T = 0. This is related to the issue of renormalization 
[111 fTl} [9]. Although we do not minimize in the trace-class in the case 
T 7^ but rather in the Hilbert-Schmidt class because the free energy is 
only coercive for the Hilbert-Schmidt norm, it turns out that the minimizer 
is trace-class nevertheless. 

The second (and related) interesting property shown in Theorem [3] below 
is that the total electrostatic potential created by the density v and the 
polarized Dirac sea decays very fast. More precisely we prove that 

Pj-v e L 1 ^ 3 ) and (p^ — v) * — G L 1 (M 3 ). 

Necessarily, the charge of p^ and the charge of the external sources have 
to be equal. More precisely the effective potential has a much faster decay 
at infinity than l/|x|, which shows that the effective potential is screened. 
In other words due to the positive temperature, the particles occupying the 
Dirac-sea have enough freedom to rearrange in such a way that the exter- 
nal sources are totally shielded. Within non-relativistic fermionic plasma 
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this effect is known as Debye- screening. Let us emphasize that in order 
to recover such a screening, it is essential to calculate the Gibbs-state in a 
self-consistent way. 

These two properties of the minimizer of the reduced theory probably 
also hold for the full BDF model with exchange term. However, like for the 
case T = 0, the generalization does not seem to be straightforward. 

The paper is organized as follows. The first section is devoted to the 
presentation of our results for the reduced model which is simpler and for 
which we can prove much more than for the general case. In the second 
section, we consider the original Hartree-Fock model with exchange term. 
We prove the existence and uniqueness of the free Hartree-Fock vacuum, 
define the BDF free energy in the presence of an external field and provide 
some interesting properties of potential minimizers. In the last section we 
provide some details of proofs which are a too lengthy to be put in the main 
text. 

Acknowledgments. M.L. acknowledges support from the ANR project 
"ACCQUAREL" of the French ministry of research. R.S. was partially sup- 
ported by U.S. NSF grant PHY-0652356. and by an A. P. Sloan fellowship. 

1 The reduced Bogoliubov-Dirac-Fock free energy 
1.1 Relative entropy 

Throughout this paper, we shall denote by & p ($j) the usual Schatten class 
of operators Q acting on a Hilbert space f) and such that tr(|<5| p ) < oo. 
The UV cut-off is implemented like in |10t [TT| [rZj [T3] in Fourier space by 
considering the Hilbert space 



with .B(0, A) denoting the ball of radius A centered at the origin. We denote 
by 'tr' the usual trace functional on &i(Sj\). Within the reduced theory, 
the free vacuum at temperature T = > is the self-adjoint operator 
acting on fy\ defined by 



Notice when T — > {j3 — > oo), we recover the usual formula |10| fTTj [9] 
/ = -D°/2\D°\. 



We assume that T > henceforth. Notice that thanks to the cut-off in 
Fourier space and the gap in the spectrum of D°, the spectrum of 7 does 




(9) 




(10) 
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not include or ±1/2. In fact, it is given by 



1 e 



-0E(A) 



1 e 
•x + 



2 i + e-W)' 2 l + e-P 



U 



-/3 1 



-/JJ5(A) 



1 



2 1 + e~P ' 2 1 + e -^( A ) 



11) 



where -E(A) = \f\ + A 2 . Also the charge density of the free vacuum 7 at 
temperature T vanishes: 



1 



2(2^) 3 



tr 



c 4 



1 



dk = 0. (12) 



>B(o,\) V l + e^°( fc ) l + e-W), 
We shall denote the class of Hilbert-Schmidt perturbations of 7 by IC: 



K := {7 € S(£ A ) I 7* = 7, < 7 < ^, 7 " 7° G 6 2 (.Sa)} • (13) 
The relative entropy reads 



#(7,7° 



tr 



^+7) (ln(i+7)-ln(i+7°)) 



+ (i- 7 ) (ln(i- 7 )-ln(i-7°)) 



(14) 



Note that since 7 G JC is a compact perturbation of 7 , we always have 
Cess(7) = fess(7°)- Hence (7(7) only contains eigenvalues of finite multiplicity 
in the neighborhood of ±1/2. Using the integral formula 



In a — In b 



a+t b+t 



dt 



00 1 1 
a+t b+t 



we easily see that Eq. (|14p is well defined as soon as 7 G IC, 7 — 7 G &i(S)a), 
since the spectrum of 7 does not contain ±1/2. 

When 7 — 7 G IC is merely Hilbert-Schmidt, we may define the relative 
entropy by the integral formula 



#(7,7° 



tr 



1 1 + 2wy' 



2 _( 7 - 7 °)i K( 7 _ 7 0)_ 1 

o 17 7J l±2u7 17 7J 1 



+ 2ii7 c 



It is clear that this provides a well defined object in IC as one has 



(16) 



for some e > 0, by (fTTj) . It is not difficult to see that (fTBT) and (fill) coincide 
when 7 — 7° G <5i(.£)a)- We shall discuss this in the appendix. But (fT6|) has 
the advantage of being well-defined for all 7 G JC, and hence we use (fT6l) for 
a definition of -ff henceforth. 
Our first result is the 

Theorem 1 (Properties of Relative Entropy). The functional 7 1— ► 
#(7,7 ) defined in $16]) is strongly continuous on AC for the topology of 
62 (^)a)- H is convex, hence weakly lower semi- continuous (wise). Moreover, 
it is coercive on K, for the Hilbert- Schmidt norm: 

V 7 G K, TH( 7 , 7°) > tr (|D°|( 7 - 7 °) 2 ) (17) 

where we recall that T = f3~ l is the temperature. 

Coercive in this context means that ^(7,7°) — > 00 if U7 — 7°|g 2 (^ A ) ~^ 
00. This follows from ( 1171 ) since \D°\ > 1. 



Proof of Theorem^ First we prove that i?(-,7 ) is strongly continuous for 
the 62 (.ft a) topology. This is indeed a consequence of the following 

Lemma 1. Let 7,7' G fC. Then we have for some constant C (depending 
on A) and all < n < 1, 

|i7(7,7 )-^(7',7 )|<^|7-7lle 2(flA) 

+ ^(!|7-7 !le 2(s , A) + ||7 / -7 || 2 62( ^ A) 

Proof. We use Formula (fTBT) and split the integrals as follows: 

1 r—l+ri rl— r) rl 

= + + ■ 

-1 J-l J-l+r) Jl-ri 




■u, 



1 + 2uy 




1 + 2n7 1 + 2wy' 1 + 2wy v ' 7 1 + 2wy' 
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and < (1 + 2wy') 1 < t] 1 as 7' G /C and —1 + 77 < u < 1 — 77. Similarly 

du — L/ 7 _ 7 o)I^R ( 0) — 1 

l + 2u7° w ' y l + 2ti7 W ' ; l + 2u7° 



<^!!7-7°||e 2( , A ) 



The other terms are treated in the same way. □ 

Convexity of 7 ^ (7, 70) is a simple consequence of the integral rep- 
resentation (]16p . In fact, the integrand is convex for any fixed u G [—1,1], 
since 

7 ^ (7 - 70) itW 7 " 70) 

= (2^ ( 2u7 " 1 " 4n7 ° + (1 + 2u7o) TT^ (1 + 2n7o) 

is clearly convex. 

Finally, we prove Formula (I17p . Consider the following function 



/(x, y) = (i + x) (ln(i + x) - ln(i + y)) + (| - x) (ln(± - x) - ln(± - y)) 

defined on (—1/2, 1/2) 2 . Minimizing over x for fixed y, one finds that 
/(x, y) > (x - y) 2 C(y) where C(y) = In f^tzfj I (2y). If we write y as 



" = Htt?-tt^J- (19) 

we obtain C{y) = /itanh(/i/2) _1 > max(|/i|, 2). Hence if y takes the form 
(fTT?|) . we deduce 

/(x,y) > max{(x - y) 2 \h\, 2(x - y) 2 } . 

Assume now that X and Y are self-adjoint operators acting on a Hilbert 
space £j, with -1/2 < X, Y < 1/2 and 

r 1/ 1 1 



2 1 1 + 1 + e- 



for some .ff. By Klein's inequality |23|, p. 330], one also has 

H(X,Y) = tvf(X,Y) >max{tr(X- Y) 2 |#|, 2tr(X- Y) 2 }. (20) 
This gives (fT71) . taking X = 7 and Y = 7 . □ 
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1.2 Existence of a minimizer and Debye screening 

Now we are able to define the reduced Bogoliubov-Dirac-Fock energy at tem- 
perature T = j3~ l . For this purpose, we introduce the Coulomb space 

C := {p G S'(R 3 ) | D(p, p) < 00} (21) 

where 

D(f,g)=4<K [ \k\- 2 7{k)9(k)dk. (22) 

We remark that the Fourier transform of Q = 7 — 7 in an L 2 -function 
with support in -6(0, A) x B(0, A). Hence Q(x,y) is a smooth kernel and 
Pq{x) = tr C 4(Q(x, x)) is a well defined function. Indeed, the map 7 G K, 1— ► 
p 7 _ 7 o G L 2 (M 3 ) is continuous for the topology of ©2(-6a)- It is easy to see 
that the Fourier transform of p 7 _ 7 o is given by the formula 

C?W = 7^4^ L,/o,^A tr c4 k^7 )(p + fc/2,p- fc/2)l dp. (23) 



7—7" W ( 27r )3/2 /|p+fc/2|<A 
V ; |p-fc/2|<A 



We also define our variational set by 

K c := {7 G /C I p 7 _ 7 o G C} . (24) 
The reduced Bogoliubov-Dirac-Fock energy reads 



^ ed ( 7 ) = TF( 7 ,7°) - aD(i/,p 7 _y.) + ^D( Pj ^o, p^ j0 ] 



(25) 



and it is well-defined on Kc by Theorem [TJ In (I25p , ^ G C is an external 
density creating an electrostatic potential — v * l/\x\. The number a > 
is the fine structure constant. The following is an easy consequence on 
Theorem [TJ 

Theorem 2 (Existence of a minimizer). Assume T > 0, a > and 
1/ £ C. T/ien J-"|? d satisfies 

V 7 G/C C , ^ ed (7) > "f D{u,v) (26) 

hence it is bounded below on Kc ■ It has a unique minimizer 7 on Kc ■ The 
operator 7 satisfies the self- consistent equation 

1/1 1 




2 V 1 + eP D i 1 + e-P°i J ' (27) 
= D° + a(p 7 _ 7 o - u) * I • 

Remark 1. When T = 0, a similar result was proved in Thm 3], but 
there might be no uniqueness in this case. 
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Remark 2. If there is no external field, v = 0, we recover that the optimal 
state is 7 — 7 = 0, and its energy is zero, by (|26l) . 

Proof of Theorem [3 Eq. (|26H is an obvious consequence of positivity of the 
relative entropy H and positive definiteness of D (•,•). The existence of a 
minimizer is obtained by noticing that is weakly lower semi-continuous 
for the topology of ©2(£>a) and C, by Theorem [TJ As is conveXi] and 
strictly convex with respect to /5 7 _ 7 o , we deduce that all the minimizers share 
the same density. Next we notice that ±1/2 ^ f (7) since the derivative of 
the relative entropy with respect to variations of an eigenvalue is infinite 
at these two points. Hence 7 does not saturate the constraint and it is 
a solution of Eq. (|2"7|) . This a fortiori proves that 7 is unique, since Z? 7 
depends only on the density /9 7 _ 7 o. □ 

Now we provide some interesting properties of any solution of Eq. (j27j) . 
thus in particular of our minimizer 7. 

Theorem 3 (Debye Screening). Assume T > 0, a > and v G C f] 

L 1 (M 3 ). Any 7 G K, that solves Eq. (|27H is a trace-class perturbation 0/7°, 
i.e., 7 — 7 G @i(^a). charge density p^_^o is an L X (R 3 ) function which 
satisfies 

I p o = I v and (p„„o -i/)* n £ LW). (28) 

This result implies that the particles arrange themselves such that the 
total effective potential (/? 7 _ 7 o — has a decay much faster than l/\x\. 

This implies that the nuclear charge of the external sources is completely 
screened. 

The proof of Theorem [3] is lengthy and is given later in Section 13.11 

2 The Bogoliubov-Dirac-Fock free energy 
2.1 Definition of the free vacuum 

When the exchange term is not neglected, the free vacuum is no longer de- 
scribed by the operator 7 introduced in the previous section. Instead it is 
another translation-invariant operator 7 that solves a self-consistent equa- 
tion. Following ideas from [T3], we define in this section 7 as the (unique) 
minimizer of the free energy per unit volume. We consider translation- 
invariant operators 7 = j(p) acting on and such that —1/2 < 7 < 1/2 
which is obviously equivalent to —1/2 < j(p) < 1/2, for a.e. p G B(0,A), in 

1 lt can indeed be proved that H(-,j°) is strictly convex but we do not need that here. 
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the sense of C 4 x C 4 hermitian matrices. The free energy per unit volume 
of such a translation-invariant operator 7 at temperature T is given by 



T t(7) = 7^( / ti c 4D°(p) 7 (p)]dp 
*B(0,A) 

« f f tr c4 [7(p)7(g)] 



rr **Wphto] ^_ TSi ) \ (29) 

JJb(o,a) 2 \p-q\ ) 



(2tt)* 

where the entropy is defined as 

S(l) = ~ [ tr C 4[(i + 7 (p))ln(i + 7(p))+(|- 7 (p))ln(|- 7 (p))]*- 
Jb(o,a) 

The free energy is defined on the convex set of matrix- valued functions, such 
that, for all p £ 5(0, A), j(p) is a hermitian 4x4 matrix, i.e. 

A:= { 1 :B(0,A)^M 4 \ 1 ( P y= 1 (p), 

- 1/2 < 7(p) < 1/2 for all p G 5(0, A)}. (30) 

Theorem 4 (The free vacuum at temperature T). For all T > and 

a// < a < 4/7T, i/ie /ree energy per unit volume Tt in $29\) has a unique 
minimizer 7 on A. It is a solution of the self- consistent equation 



7 



_ 1 ( 1 1 



n- n - no _ n . 7°(s,y) 



2Vl + e /3 f;\ l + e _/ "W (31) 

-D^o - 

Furthermore, 7 /ias t/ie /orm 



7°(p) = /l(H)«-p + /o(bl)^ (32) 

with /o,/i < a. e. on 5(0, A) and D^o satisfies 

\Dy>\ > \D°\. (33) 

Here and in the following, we shall identify operators with their integral 
kernels for simplicity of the notation. That is, the last term in the second 
line of (]3ip denotes the operator with integral kernel given by ^ 7 , 
where (7—7°) (x, y) is the integral kernel of the translation- invariant operator 
7 — 7 (it is a function of a; — y). 

Remark 3. The assumption a < 4/tt guarantees that the functional (|29h is 
bounded from below, independently of the UV cutoff A, which is arbitrary 
in this paper. This is a consequence of Kato's inequality. For a > A/tt this 
is not the case [T7] . 

For comparison, we note that in the non-interacting case a = 0, the 
functions /i(|p|) and /o(b|) appearing in Theorem 2] are given by 

/i(bl) = /o(bl) ' / ' 1 
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A similar result was proved in the zero temperature case in [T5]. As in 
j, it is possible to justify the introduction of Tt by a thermodynamic limit 
procedure. The proof of Theorem 0] is given in Section 13.31 
Like for the reduced case, we have that 

a(7°) C 

for some e > 0. This can be seen from (|3"3"|) and the fact that D^o is a bounded 
operator on Sja due to the presence of the ultraviolet cut-off. Notice also 
that we have formally p^o = by ([32]) . as in (fT2~|) , 

2.2 The external field case 

As in Section 11.2} one can consider the Bogoliubov-Dirac-Fock energy with 
an external field. It is formally obtained by subtracting the infinite free 
energy of the free vacuum at temperature T > from the free energy of our 
state 7. This procedure can be justified like in [14J by a thermodynamic 
limit procedure. Using the same notation as in Section 11.21 the Bogoliubov- 
Dirac-Fock free energy reads 

.F T (7) = TH(j,j°) - <xD(z/,p 7 _ 7 o) + ^£>(p 7 _ 7 o,p 7 _ 7 o) 

a ff tlC ^-^ y )ld xd y, (34) 



2 J J \x — y\ 

where H is the relative entropy defined like in Section 11.11 Like for the 
reduced case, we see that the functional Ft is well-defined on the following 
convex set 

ICc := {7 G B(fi A ) I 7* = 7, ~\ < 7 < \, 7 - 7° G © 2 (^a), P 7 -f» G C 

(35) 

Note that although the function 7 1— > ^(7,7°) is convex, Ft is not a convex 
functional because of the presence of the exchange term. This is of course a 
great obstacle in proving the existence of a minimizer, and we have to leave 
this as an open problem. Following the method of Theorem [TJ we shall show 
that 

V 7 G/C C , TF( 7 ,7°) >tr(|D 7 o|( 7 -7°) 2 )- (36) 
With the aid of this inequality we can prove the 

Theorem 5 (Minimizer in External Field). Assume that < a < A/tt 
and that T > 0. We have 

V 7 g£ c , Fr(i) > ~D(v,v) (37) 
and hence Tt is bounded below on ICc- 
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Assume that 7 G /Cc is a minimizer of Tt- Then it satisfies the self- 
consistent equation 



1 



1 



7 



1 



2 I 1 + eP°i 1 + e~P D i J ' 



Z? 7 := DfM + a(p 7 _ 7 o - z/) * 



-1 _ „,(7-7°)fa«) 



F~ 1/1 



TOf/i Z? 7 o defined in [31]) . It is unique when 



(38) 



°^l{ 1 - a (f/^ + ' 1/62lI/6 ) D (-) 1/2 } ^ < 39 > 
The proof of Theorem is provided in Section 13.41 



3 Proofs 

3.1 Proof of Theorem H 

Let 7 be a solution of 

1 



7=o 



1 



2 V 1 + 1 + e-^ D n 

Z? 7 := D° + a(p 7 „ 7 o —u)*\- 



(40) 



For the sake of simplicity, we define p := /? 7 _ 7 o — and V = a(/3 7 _ 7 o — v) * 
I • I" 1 . Note that W G L 2 (R 3 ) as p G C, hence V G L 6 (M 3 ). Following [TH 
p. 4495], we may use the Kato-Seiler-Simon inequality (see [20J and [2T| 
Thm. 4.1]) 

Vp > 2, |/(-iV) 5 (x)|| ep(L2(R3)) < (2vr)^ 3 /f |M| LP(R3) l/l^a) (41) 
to obtain 



V 



\D°\ 



< 



V 



\D°\ 



< a 



S 6 (fi A ) 



<C\ 



lie • 



This shows that \D^\ < (1 + ctC \p\c)\D Q \. Thanks to the cut-off in Fourier 
space, we deduce that D 7 is a bounded operator or Sja- Recall Duhamel's 
formula 

e PD, = e (lD» + p f 1 e tpD,y e (l-t)l3D° dt (42) 

Jo 

Denoting K := (3 Jq 1 e t0D y y e (i-t)pD° ^ and us i ng ([^j), we have 
K = K + K' 

:=P f l dte^ D °Ve^ D ° +I3 2 f dt f ds e s ^Ve^ D °Ve^ D ° . 
Jo Jo Jo 
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We obtain for the self-consistent solution 

7-7' 



.,0 _ 1 



1 + 1 + e pn° 

1 K 1 1 K , 1 

1 + e /3D0 ^0 x + g/3D o 1 + e/3D o * x + e/3D0 

+ l + e /3D<» X 1 + e /3D 7 K 1 + e/ 3D0 ( 43 ) 

which we write as 7 — 7 = A + B where 

1 1 yi e^o° e {i-t)f3D° 



As 7 G L b (R d ) and L> 7 is bounded, using the cut-off in Fourier space and the 
Kato-Seiler-Simon inequality (|4ip. we have K G ©6(^a)- Hence we obtain 
that if' G 6 3 (£a) and B G 6 3 (£a)- 

The next step is to compute the density of A. The kernel of A is given 

by 

pi tpD°(p) _ Jl-t)f3D°(q) 



Using (|23|) . we obtain 



aC(\k\) 



paw = 



where 

op pi r tl3D°(p+k/2) Jl-t)0D°(p-k/2) 

|p-fc/2|<A U 



1 + e /3D°(p+k/2) 1 + e f3D°{p-k/2) 



Inserting this into the self-consistent equation (|43p gives 



(44) 



m = -m- ^p-m +PBik) (45) 

or, equivalently, 

p(k) = bi(k)(-v(k)+p£(k)), (46) 

and 

V(k) = 4Trb 2 (k)(-v(k) + p%(k)), (47) 

where 

fe l : = f i,i9 ^^/^u l and & 2 := Cttio ^TTTiT ^ > ( 48 ) 

Vlfc| 2 + «^(l^l)y vl^l 2 + «c(|fc|)y 

with F~ l denoting the inverse Fourier transform. Our main tool will be the 
following 
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Proposition 1 (Properties of 61,62)- The two functions bi(x) and 62(2), 
defined in (gSJ and ([g]). 6e/on# to L^IR 3 ). 

We postpone the proof of Proposition Q] to Section EO and first complete 
the proof of Theorem First we claim that ps £ L 3 (R 3 ). To see this, we 
take a function £ £ L 3 / 2 (M 3 ) n C^QR 3 ) and compute 

\tr(BO\ =|tr(51 B(0iA )b)as(o,A)(p))l 

^ I^IssWa) ll 1 B(O 1 A)(p)ClB(0 I A)(p)| e3/2(i5A) • 

Writing £ = |f| 1/2 sgn(f)|£| 1/2 and using the Kato-Seiler-Simon inequality 
dUD twice in &s(Sja), we obtain | tr(B£)\ < C ||-B|le 3 (£ A ) ||£|| L 3/2( M 3) where C 
depends on the cut-off A. This proves by duality that ps £ L 3 (M 3 ). 

Next we use a boot-strap argument. As v £ L 1 (1R 3 ) and ps £ L 3 (1R 3 ), 
we get from (|47p and Proposition [1] that V £ L 3 (M 3 ). Inserting in the 
definition of K' and using ([4"T]) once more, we obtain that K' £ ©3/2 C^a)) 
hence B £ 6 2 (K 3 ) and p_e £ L 2 (M 3 ). Using again (|37j) and Proposition [Q 
we get that V £ L 2 (R 3 ), hence 5 £ &i(S)a) and p B £ L^R 3 ). This finishes 
the proof of Theorem El by (|46p . (I47D and Proposition [TJ □ 



3.2 Proof of Proposition [T] 



The proof proceeds along the same lines as in the Appendix of [9]. In the 
following we shall denote by Pq and Pq~ the projection onto the positive and 
negative spectral subspace of D°, respectively. As multiplication operators 
in momentum space, 



1 + 



a ■ p + P 

W) 



1 



a -p + /3 



The function C in (I44p can be written as 



p+fc/2|<A 
p-fc/2|<A 



(fox 



(it tr C 4 



tr C 4 



e t/3E(p+k/2) e (l-t)0E(p-k/2) 
1 + e PE(p+k/2) l + e f3E(p-k/2) 

e t/3E(p+k/2) e -(l-t)PE(p-k/2) 



P+(p + k/2)P+(p-k/2) 



1 + e /3£(p+fe/2) 1 + e -f3E(p-k/2) 



P+(p + k/2)P^(p-k/2) 



(49) 
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Hence 



C(\k\) 



(3E(p+k/2) _ J3E(p-k/2) 



1 



1 

^2 J\ P+ k/2\<A i + e PE(p+k/2) E {p + k/2)- E(p - k/2) 1 + e^(P- fc / 2 ) 

\p— fc/2|<A 



1 + 



(p + k/2) ■ (p - k/2) + 1 
E(p + k/2)E(p-k/2) 



dp 



+ 



e 0E{p+k/2) _ e -/3E(p-k/2) 



7T 2 J\P+k/2\<k 1 + e 0E{P+k/2) E ( p + k / 2 ) + £(p _ 1 + e ~fiE(p-k/2) 

\p—k/2\<A 



{p + k/2) • (p - fc/2) + 1 
E(p + k/2)E(p-k/2) 



dp. 



For the sake of clarity, we denote by Ci(|/c|) (resp. C 2 ( | A; | ) ) the first (resp. 
second) integral of the previous formula. By the monotonicity of the expo- 
nential function, it is easily seen that Ci(|fc|) > and C2 ( I ^ I ) — 0- 

The next step is to simplify the above integral formula. We follow a 
method of Pauli and Rose [19] which was recently used in the appendix of 
[9]. After two changes of variables, we end up with 



dpi) 



7r|A;| 



(3w(k,z) S i n h(/3u) 



dz / dv- 

l l + e P(w(k,z)+v) 



1 



X 



z w(k, z) 1 



7T 



^-z A {\k\) 



1 + e p(w(k,z)-v) (l_ z 2)3 

e /3(B(A)- z ) S mh(/3t;) 1 



dz dv 

l 14. e/9(£(A)-z-N)) 



(W) - .) 2 - M-) , 



1 + e/9(B(A)-*-«) 

(50) 



ir\k\J J 1 + e WM+») 1 + ffi (1 _ z 2) 

|fc| 2 _ 



1 + e p(v-w(k,z)) ^4 ; 1 _ z 2 

8 r l 4 z ^\ k \) r* e 0v sinh(/?(£(A) - z)) 

+ W\'o Jo £(A) - i 



1 



1 + e /3(v+z-E(A)) y 4 
In the above formulas we have used the notation (as in |9j) 



Za(t) 



Vl + A 2 - y 7 ! + (A - r) 2 



(51) 
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Note that is a decreasing C°° function on [0, 2A] satisfying Z\(0) 
A/E(A), Z\(2A) = 0. We have also used the shorthand notation 



"(M)=y i+i Ty )/4 . 

All integrands of the above formulas are real analytic functions of r = \k\ on 
a neighborhood of [0, 2A]. Also all the integrals vanish at k = 0. We deduce 
that C\ and C2 are smooth functions on [0, 2A]. Using Za(2A) = 0, one also 
sees that Ci(2A) = C((2A) = C 2 (2A) = C 2 (2A) = 0. A Taylor expansion of 
the first integral of C\ yields 

4 r E{ - k) t 2 dt 

d(0) = -P / m > 0. 

The end of the proof of Proposition Q] is then the same as in [9j Prop. 17]. 
First we notice that as C(r) is bounded and has a compact support, b\ and 
62 are in L°°(M 3 ). We now prove that they decay at least like |x| -4 at infinity 
meaning that they also belong to L X (R 3 ). To this end we write for b = b% 
or = 62 the inverse Fourier transform in radial coordinates: 

1 f 2A - 

VxGM 3 \{0}, b(x) = (rb(r))ean(r\x\)dr. (52) 

V27r|x| Jo 







Integrating by parts and using 6(2 A) = b'(2A) = yields 

Vx £ R 3 \ {0}, b(x) = * - I 2A6"(2A) cos(2A|a;|) - 26' (0) 

V2vr|x| 4 \ 

2A _ \ 

(rb)^\r) cos (r I a; |) dr J . (53) 

This completes the proof of Proposition [U □ 
3.3 Proof of Theorem H 

The proof is inspired by ideas from [II]. We denote I := inf 7e ^4 7^(7). We 
start by introducing the following auxiliary minimization problem 

J=infT T ( 7 ) (54) 

7S£5 

where B C A is given by 

B:={ 7 e A, 7 (p) = /i(b|)a • p + f (\p\)P, fo, A < 0} . (55) 
Lemma 2. There exists a minimizer 7 6 B for (I54p . 
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Proof of Lemma [B The functional Tt is weakly lower semi-continuous for 
the weak-* topology of L°°(B(0, A)). This is because —S is convex and the 
exchange term is continuous for the weak topology of L 2 (B(0, A)) as shown 
in [13]. Also B is a bounded closed convex subset of L°°(B(0, A)). Hence 
there exists a minimum. □ 



Lemma 3. Let ^° £ B be a minimizer of (I54p . Then there exists an e > 
suc/i tfwti | 7 °| < 1/2 - e. 

Proo/. For x £ [1/2,1/2], 

s(x) := (| + x) In (1 + x) + (| - x) In (± - x) 

is an even function of x. Because of the special form of 7 , we have 7 °(p) 2 = 
|| 7 (p)|| 2 -fc 4 for all p E -B(0, A), where || • || denotes the matrix norm. Hence 



V 7 e£, 5( 7 ) = -4/ s (|| 7 (p)||)dp. 

JB(0,A) 

The derivative of s is infinite at x = 1/2 and the derivative of the terms of the 
first line of (|29p stays bounded. It is therefore clear that {p G B(0, A) | 1/2 — 
e < 1 1 7 (p) 1 1 < 1/2} has zero measure for e small enough. □ 

Let us now write the first order condition satisfied by 7 . Since ||7°(p) || < 
1/2 — e for some e small enough, we can consider a perturbation of the form 

7 (p) = f( p ) + f ( 5l (|p|)a • p + ffo(|p|)/3) 

with <70)9i < and i > small enough. We obtain 

r,0/ 



/ tr C 4 
JB(0,A) 



^ 0(P) + Tln 1/2 -70S ) (5l(H)Q ' P + 5 ° (H)/?) 



for all gi,g < 0. 

We notice that the function x 1— > In f 1/2-3; ) * s °dd, hence 



dp > 
(56) 



with sgn(7) = 7/I7I. We obtain that 



ln( ^i^li) = ^ )F( ^ w||) 



where F(x) = In ^ j^-x ) Z 3 "' ^ n ^ e °ther hand, we can write 

=d 1 {\p\)a-p + d Qp\)/3 
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where d\ and do are given by [HI Eq. (72)-(73)]. Using /i,/o < 0, we 
immediately see that 

di(H)>l and do(|p|)>l, (58) 
which in particular proves that 

[|^o(^)[| > [|J^°0)][ > |^|. (59) 

All this gives 

D^+T\n^^ = {d 1 (\p\)+Tf 1 (\p\)F(\\f(p)\\)) a -p 

+ (do(N) + T/ (H)F(||7°(p)||))/3. (60) 
Inserting this in (|56p . we obtain the first order conditions 



diw+Tfi(\pm\\^(p)\\)<o 

do(bl)+T/o(b|)F(||7°(p)||)<0. 

In particular, because of ([581) we infer that 

/i(N)F(||7°(p)||) < -1/T 
M\P\M\\T°(P)\\) < -i/T. 



(61) 



(62) 



As -F(||7 (p)||) > and /o > — ||7°(p)||, we obtain from (|62l) the inequality 
l7(p)|F(ll7(p)ll) > Hence 

This inequality means that /o and f\ cannot vanish simultaneously. But we 
can indeed prove that each of them cannot vanish, as expressed in the 

Lemma 4. Let 7°(p) = /i(|p|)a -p + fo(\p\)P be a minimizer of (|54j) . Then 
there exists an e > such that 

fo < -e and fx < -e. 

Proof. By Lemma [3] we know that [I7 (p)\\ < 1/2 — e for some e > 0. By 
(|62l) and the monotonicity of F we obtain /& < — T pHT^e) for = 0, 1. □ 

Lemma 5. Let j° £ B be a minimizer of (|54|) . TTien solves the self- 
consistent equation 

7 ° = 2 (l + e ^o " 1 + e -^o) ' (64) 
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Proof. As the constraints are not saturated by Lemmas [3] and HJ we obtain 
that the derivative vanishes, i.e. 



d l (\p\)+Tf 1 (\p\)F(\\f(p)\\) = 
do(bl)+T/o(b|)F(||7°(p)||) = 0. 



(65) 



which means that 

1/2 + 7° 
Z^o+Tln 1^0=0- 



Hence 7 solves (|M|) . □ 

Now we prove that the operator 7 defined in the previous step is the 
unique minimizer of Tt on the full space A defined in (|30l) . not merely on 
the subset B in (fo3|) . We have 

r T(7 ,-r TR ») = th^^SL. ff ^IO-f)W(7--WW) ^ 



(2tt) 5 JJb(o,A) 2 \P~1 
where H is the relative entropy per unit volume 

iJ( 7 , 7°) = (2vr)- 3 / tr C 4 [ (i + 7) (In (± + 7) - In (± + 7°)) 

+ (i-7)(ln(i-7)-ln(i-7 ))]^- (66) 

We shall use the important 

Lemma 6. For H in 166\) the inequality 



Ttf (7,7°) > (2vr)" 3 / tr C 4 |D*o(p)|( 7 (p) - 7 °(p)) 2 ^ (67) 



B(0,A) 



holds for all 7 € .4.. 



Proof. This is a simple application of (|2Up . taking X = ^y(p), Y = 7°(p) and 
integrating over the ball B(0, A). □ 

Using Lemma [6] and the formula 



(2vr)2 JJ mA)2 \p-q\i " 2 



where j(x) is the Fourier inverse of the function 7(p), we find 



rr(7)-Tr(7°)> (2i)- s ( / tr c , |D,,(j>)|(7(f>) - 7°(p)) 2 * 

W ' A> _W t r c ,|(7-7°)M|^ m 
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for all 7 £ i. We now use ideas of [TJ QUI HI]- Kato's inequality |x| 1 < 
7r/2|V| gives 

« /• t.- c ,|(7 - 7°)Wl^ ^mr f _ 
2 Jr3 Fl 4 Jb(o,a) 

By (f5"9"l) we deduce 

T T ( 7 ) - T T (7°) > (1 - vra/4)(2vr)- 3 / tr C 4 |Z) fJ (p)|( 7 (p) - f(p)) 2 dp. 

Jb(o,a) 

Hence 7 ° is the unique minimizer of Tt on A when < a < 4/tt. This 
completes the proof of Theorem (5J □ 

3.4 Proof of Theorem [5] 

The lower bound (I37p is obtained by following an argument of [U [10]. By 
(|36|) we have for all 7 G /Cc* 

Jt<7) > tr (|^|( 7 - 7°)*) - I// trc,| ^;_y-'' )l ' to» 

+ 2 D {P~,-^ ~ P7-70 -v)- 2 D i u > v )- ( 69 ) 
By ([S3]) together with Kato's inequality |x| _1 < (vr/2)|V| < (vr/2)|Z?°|, 
tr C 4|(7-7°)(^y)| 2 , , / tt 



f - y| 



dx^<-tr(|^o|(7-7 ) 2 ) 



which yields (|37p when < a < 4/vr. 

Assume now that 7 is a minimizer of Tt on ACc- The proof that it satisfies 
the self-consistent equation (|38H is the same as in the case of the reduced 
BDF functional in Section 11.21 Note that because of (|69[) and inf^c — ^' 
we obtain 

2 a\ ff tr C 4 K7 — 7 ) fx, y)| 2 , , a , ^ a' , s 

n~2j JJ \x-y\ d V + 2 D{P ^°> P ^ 0) ~ 2° M - 

It was proved in [11, p. 4495] that this implies (under the condition (|39|) ) 
that 

\D~\ > d- 1 \D°\ (70) 



'7 
with 



d = I 1 — a I — 4 
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Next we fix some 7' G JCc and use that 
H{ 1 ',f) = H{ 1 ,f) + H( 1 \ 1 ) 



+ tr 



(V- 7 ) In I 



1 + 7 V 



In 



7 



1-7°. 



(71) 



Inserting Eq. (I38p for our minimizer 7 and Eq. (|31f) for 7 , we obtain for 
any 7' G ICc the formula 



tr C 4 K7 - 7 ) (x, 



F - y\ 



-dxdy. 



(72) 



We may use one more time (|15j) and the self-consistent equation (|38j) to 
obtain 

TF(7',7)>tr(^ 7 |(7'-7) 2 ). 
By (jTUj) and Kato's inequality as before, we eventually get 

tr C 4|(7-7°)(z,y)| 2 



a (2a 
Trii) > JF T ( 7 )+-I}(py_ 7 ,p Y _ 7 )+( — - 2 



F - y| 



dxdy. 



Hence we obtain that any minimizer is unique when and/4 < 1, as stated. 
Let us remark that the expression in last term of (|7ip is indeed a trace-class 
operator. It would, however, have been sufficient to choose 7' as trace class 
perturbation of 7 and conclude the rest by a density argument. □ 



A Appendix: Integral Representation of Relative 
Entropy 

Here we shall prove the claim made in Section fTTTI that (I14p and (|16p coincide 
as long as 7 — 7 G and the spectrum of 7 does not contain ±1/2. 

From the integral representation (|15|) . we have 



a(lna — ln6)= / (a — b) dt . 

a+t b+t 



We split the first factor as 



a b t , , , 1 

+ ir—M- h > 



a + t b + t b + t K 'a + t 
and obtain 



o(ln a — In b) = I 
Jo 



00 t 1 1 

(a - b) (a - b)- dt 

b + V 'a + V 'b + t 

00 b , ,,1 



+ / - 7 (a-b)- dt. 

/o b + V 'b + t 
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Now if a — b is trace class and the spectrum of b is contained in (0, oo), then 

r°° b 1 f 00 b 
tr / (a - b)- dt = tr / (a - b)— ^dt = tr(a - b) . 



6 + t b + t J (b + t) 

If we apply this reasoning to a = 1/2 + 7, b = 1/2 + 7 and to a = 1/2 — 7 
and 6 = I/2 — 7 , respectively, we thus obtain 

(r°° t 1 1 

/ T T7T7 (7-/) ?T ^(7-/)r T ^ 

Changing variables from t to it = 1/(1 + 2t) in the first integral and u = 
— 1/(1 + 2t) in the second integral, respectively, we arrive at the integral 
representation (|16p . 
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